Quantum phase transition and excitations of the Tavis-Cummings lattice model 
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The enormous progress in controlling quantum optical and atomic systems has prompted ideas 
for new experimental realizations of strongly correlated many-body systems operating with light. 
These systems consist of photons confined in optical cavities, which interact strongly with atoms or 
atomiclike structures. Due to the interaction between the two particle species optical nonlinearities 
appear, leading to a quantum phase transition from Mott to superfluid phase. Here, we address the 
Tavis-Cummings lattice model, which describes light-matter systems containing multiple atomiclike 
structures in each cavity. In particular, we investigate the phase boundary delimiting Mott from 
superfluid phase and the elementary excitations of the two-dimensional Tavis-Cummings lattice 
model in dependence of the number of atomiclike structures per cavity. In order to obtain the 
results we employ the variational cluster approach. We evaluate spectral functions and densities 
of states of both particle species, which allows us to characterize the fundamental excitations of 
light-matter systems. These excitations are termed polaritons and are superpositions of photons 
and atomic excitations. We introduce polariton quasiparticles as appropriate linear combinations 
of both particle species and analyze the weights of their constituents. Our results demonstrate 
the dependence of the quantum phase transition and the elementary excitations on the number of 
atomiclike structures per cavity and provide thus valuable insight into the physics of light-matter 
systems. 

PACS numbers: 71.36.-|-c,42.50.Ct,67.85.Dc,64.70.-p 



I. INTRODUCTION 

The push towards the experimental realization of 
quantum computers lead to incredible advances in the 
fields of quantum optics and atomic physics. Unprece- 
dented experimental control in these fields allowed to 
envision new realizations of strongly correlated many- 
body systems, which operate with light.—"— Confined light 
modes in coupled cavity arrays arc able to tunnel between 
adjacent sites and thus propagate on a lattice of cavi- 
ties. Strong correlations in turn can be observed when 
a repulsive interaction between photons is present.— This 
repulsion, which is termed optical nonlinearity, can be 
achieved by coupling the light modes to matter in the 
form of atoms or atomiclike structures present within 
each cavity. The interaction between the light modes 
and atomic like structures is achieved by means of dipolc 
coupling. In theory there exist two major schemes to 
obtain this interaction. The first is to model the atom- 
iclike structures by two-level systems, leading to an in- 
teraction of the Jaynes-Cummings type^"^ whereas the 
second approach is based on electromagnetically induced 
transparency^ and uses four-level systemsi^ In both cases 
optical nonlinearities between photons arise as the en- 
ergy for adding two photons to the cavity is larger than 
twice the energy needed to add one photon. This behav- 
ior leads to intriguing experiments such as the photon 
blockade effect)^ where only one elementary excitation is 
present in the cavity at the same time. The elementary 
excitations in light matter systems are termed polaritons. 
Polaritons are superpositions of both particle species, 
namely photons as well as excitations of the atomiclike 
structures. Following these considerations and arranging 



multiple cavities on a lattice leads to a strongly correlated 
phase in which photons arc involved. As a result, light- 
matter systems exhibit a quantum phase transition from 
Mott phase where polaritons are localized in the cavities 
to superfluid phase where polaritons are delocalized on 
the whole latticeji"— 

Up to now an experimental realization of light-matter 
systems is still missing, however, there are several promis- 
ing approaches such as quantum dots grown in photonic 
crystal cavities, transmission line cavities and toroidal or 
disk shaped cavities i^ The advantage of light-matter sys- 
tems is that they are of mesoscopic size and thus allow for 
a direct addressability of each lattice site and good ex- 
perimental control on the system parameters. Exhibiting 
these valuable properties light-matter systems might be 
used as quantum simulators for other strongly correlated 
many-body systems such as the Bose-Hubbard model^ 
or find their applications in quantum information pro- 
cessing. 

Experimentally it might be more feasible to place mul- 
tiple atomiclike structures within one cavity. Therefore it 
is important to study light-matter systems which contain 
more than one two-level system per cavity. The theoreti- 
cal model describing a single cavity with N two-level sys- 
tems is termed Tavis-Cummings modelJ^iii In the case 
of A^ = 1 it reduces to the Jaynes-Cummings modeli^ 
Light-matter systems with coupled Jaynes-Cummings 
cavities have already been investigated to some detail in 
Refs. [2 ll2| - [2(]| and will thus not be addressed here any- 
more. However, systems of coupled-cavity arrays with 
more than one two-level system per cavity have been 
rarely studied in literature. In particular, the quantum 
phase transition has been investigated on mean field level 
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FIG. 1: (Color online) Illustration of the TCL model. Blue 
bubbles represent two-level systems and red wavy arrows pho- 
tons. Photons and two-level system interact via dipole cou- 
pling. 



by N. Na et al. in Ref.l2ll. and in one-dimension with den- 
sity matrix renormalization group (DMRG) by D. Rossini 
et al. in Ref. [ij. In the present paper, we investi- 
gate the quantum phase transition and the elementary 
excitations — the polaritons — of coupled Tavis-Cummings 
cavities arranged on a two-dimensional lattice. In partic- 
ular, we evaluate the phase boundary delimiting Mott 
phase from superfluid phase for different number of two- 
level systems per cavity. Furthermore, we study spectral 
properties of both photons as well as atomic excitations 
which in turn allows us to characterize the polaritonic 
properties of the system. In order to evaluate the quan- 
tum phase transition and the spectral excitations we em- 
ploy the variational cluster approach)^ 

This paper is organized as follows. The Tavis- 
Cummings lattice model is introduced in Sec. |IT1 De- 
tails on the variational cluster approach are covered in 
Sec, mil In Scc.|IV]we discuss the results obtained for the 
phase boundary delimiting Mott from superfluid phase. 
Section |V] is devoted to the excitations of the Tavis- 
Cummings lattice model. Here, we present spectral func- 
tions and densities of states of both particle species and 
discuss polariton quasiparticle excitations. Finally, we 
summarize and conclude our findings in Sec. I VII 



II. THE TAVIS-CUMMINGS LATTICE MODEL 

A single cavity at lattice site i containing A^ two- 
level systems is modeled by the Tavis-Cummings (TC) 
Hamiltonianr^iii 

Hl^' = CO, a] a^ + e (^f + ^ ) + .g (a^ S+ + c4 Sr) , (1) 

where uJc is the resonance frequency of the cavity, e is 
the energy spacing of the two-level systems, and g is 
the atom-field coupling constant (see a single cavity in 
Fig. [1] for illustration). The operators a] and a^, respec- 
tively, create and annihilate photons in the cavity i. The 
ensemble of two level systems is described by collective 
spin operators S*" — X)i/=i "'"j' where a G {z, +, — }, 
and a^j = crj^ ■ ± ia^ are the spin raising and lower- 
ing operators. When starting from the dipole interac- 
tion between photons and two-level systems two addi- 



tional terms arise in the Hamiltonian, which are propor- 
tional to a\Sf and a^ S~ . However, for the condition 
\^c — e| <C oJc, e these terms are fast oscillating in com- 
parison to a^ S^ and a\ S~ and can thus be neglected, 
which is known as rotating wave approximatiouj^^ The 
difference between the resonance frequency of the cavity 
UJc and the energy spacing of the two-level system e is 
termed detuning A = Wc — e. As a consequence of the 
rotating wave approximation the total number of excita- 
tions fii = a\a^ + Sf + N/2 is conserved. Additionally, 
the total spin S"^ is a conserved quantity as well. The 
ground state of the TC model is always in the sector of 
maximum spin S = N /2^ which will thus be considered 
in further calculations. 

The full model consists of M coupled cavities, which 
are arranged on a lattice. Therefore we refer to this 
model as Tavis-Cummings lattice (TCL) model. The 
TCL Hamiltonian is given by 



H 



TCL 



-tj2aU,+J2^^''-f'^p 



(2) 
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where the first term allows photons to tunnel between 
cavities i and j. The tunneling strength t is given by the 
overlap integral of the photonic wave functions, which is 
considered to be nonzero only for nearest-neighbor sites 
i and j. The restriction to nearest neighbors is denoted 
by the angle brackets (■ • • ) around the summation in- 
dices. The second term of Eq. ([2]) describes the physics 
of the individual cavities and the last term controls the 
average particle number of the system, where fj. is the 
chemical potential and Np = ^^ hi is the total particle 
number. Figure[l] illustrates the TCL model. For the 
TCL model the total particle number Np is conserved as 
well as the total spin S"^ of each cavity. As in the case 
of the Jaynes-Cummings lattice model^ the TCL model 
can be rewritten as 



(h j) i 

+ 3 E(°^ ^^ + «I ^r) - (a^ - ^c) Np . (3) 



In the forthcoming discussions and calculations we use 
the dipole coupling g as unit of energy. Therefore the 
physics of the TCL model depends only on three inde- 
pendent parameters, namely, the hopping strength t, the 
detuning A, and the modified chemical potential /x — uJc- 



III. THE VARIATIONAL CLUSTER 
APPROACH 

In order to evaluate the boundary of the quantum 
phase transition from Mott phase to superfluid phase 
and to investigate the excitations of the TCL model we 
employ the variational cluster approac h^^i^^" — (VCA), 
which yields the single-particle Green's function G(k, uj) 



of the physical system H^'~^^ . VGA has been previously 
applied to light-matter systems in Refs. llSlilil . 

The basic idea of VGA is that the grand potential 17 
can be expressed as a functional of the self-energy Z and 
that Dyson's equation for the Green's function is fulfilled 
at the stationary point of ri[Z]. To be able to evaluate 
ri[II], the unknown self-energy Z of the physical system 
is approximated by the self-energy of an exactly solv- 
able, so-called reference, system. In VGA the reference 
system is chosen to be a cluster decomposition of the 
physical system, which means that the system of size M 
is divided into clusters of size L with, however, different 
single-particle parameters x as compared to the physical 
system. Due to the approximation in the self-energy the 
functional Vl\T.] becomes a function of the single-particle 
parameters x of the reference system 

r!(x) = f7'(x) + Tr ln(-G'(x)) + Tr H-{G^^ - Z(x))) , 

(4) 
where quantities with prime correspond to the reference 
system and Gq is the noninteracting Green's function. 
The stationary condition on 51 (x) now reads 



dn{y) 
dx 







(5) 



which can be evaluated numerically by varying some or 
all single-particle parameters of the reference system. In 
Refs. J4, 26. 27. it was pointed out that the particle num- 
ber of a certain particle species is thermodynamically 
consistent only if the corresponding chemical potential 
is considered as variational parameter. However, there 
is a formal difficulty when considering, for instance, the 
chemical potential ji or the energy spacing of the two- 
level systems e as variational parameter since both cou- 
ple to atomic excitations, which cannot be regarded as 
noninteracting particles. Generally a variation in these 
terms is not allowed within VGA, yet, this subtlety can 
be circumvented by mapping the atomic excitations onto 
hard-core bosons. 



h\ 
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\-S), -> |0).,, 1-^+1), ^ |1)., ...|5), ^ |iV), . 
With this mapping the TGL Hamiltonian is given by 
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TCL 
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N 



L-^E 11(^1^.--) 



U^oc 2 



i u=0 



(6) 



where /± - p{S + 1) - {b\ 6, - S){b\ h^ - S ± 1)). In 
Eq. ([5]) wc have formally added the hard-core constraint 
by introducing the infinite interaction U . Due to the 



fact that the Hamiltonian consists now of noninteracting 
particles the variation in the chemical potential /x, which 
is needed to guarantee that the total particle number Np 
is thermodynamically consistent, becomes possible. 

In our calculations we varied the on-site energies and 
the hopping strength. In fact only two of the three on-site 
energies Wc, e, and /x are linear independent and therefore 
a combination of any two of them used as variational 
parameters yields identical results. In particular, we use 
x = {u>c, e, t} as variational parameters, which ensures 
thermodynamic consistency for the particle number of 
both particle species and thus consistency for the total 
particle number as well. 

The reference system defined on a cluster of size 
L is solved using the band Lanczos method j^^'^^ The 
initial vector of the iterative band Lanczos method, 
which is used to evaluate the particle term of the 
Green's function, consists of 2L elements and is given by 

{a\ l^o) , a\ l^o) . . . ai |Vo) , h\ IV-o) , 4 IV'o) ■ • ■ &I IV-o)}, 
where jfAo) is the A'p-particle ground state. In order to 
evaluate the hole term of the Green's function the cre- 
ation operators are replaced by annihilation operators. 
Using the bosonic Q-matrix formalisni^S we evaluate 
the grand potential and the single-particle Green's 
function of the physical system. As the TGL model 
consists of two distinct particle species, we extract a 
Green's function for photons G^''(k, w) = G i{uj) 

k k 

and a Green's function for the atomic excitations 
^■"^^(k, LLj) = G, ^t(w). From the Green's functions we 

k k 

evaluate the single-particle spectral function 



1_ 



A^(k, w) = ImG'Xk, w) , 

IT 

the density of states 

iV-H^ /"A-(k,c.)dk^l^A-(k,c.) 

•' k 

and the momentum distribution 



n^(k) 



A^(k, uj)duj , 



(7) 



(8) 



(9) 



where x can be either ph for photons or ex for atomic 
excitations. 



IV. QUANTUM PHASE TRANSITION 

The nonlincarities which arise due to the coupling of 
the photons to the ensemble of two-level systems lead to 
a quantum phase transition from Mott phase to super- 
fluid phase. The elementary excitations in light-matter 
systems — the polaritons — are linear combinations of pho- 
tons and atomic excitations. The Mott phase is charac- 
terized by integer polariton density, zero compressibility, 
and a gap in the spectral function. Intriguingly, polari- 
tons in Mott phase are localized in cavities, which in turn 
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FIG. 2: (Color online) Phase boundary of the two-dimensional 
TCL model evaluated for zero detuning A = 0, reference sys- 
tems of size L, and A'^ two- level systems per cavity, where in 
(a) N = 2, (b) N = 3, (c) iV = 5, and (d) N = 10. 



TABLE I: Ratio wi /wi of the width of Mott lobe 1 with po- 
lariton density Up = 1 and Mott lobe i with polariton density 
Up = i for A'^ two-level systems per cavity. 



N 



'Wl/W2 Wl/Wi 



2 
3 
5 
10 



1.18 
1.05 
1.02 
1.00 



2.84 
1.38 
1.09 
1.02 



Wl/Wi 



4.83 
2.26 
1.25 
1.04 




FIG. 3: (Color online) Critical hopping strength ratio 
t*{N)/t* (2) for the first Mott lobe in dependence on the 
dimension D. Results for one dimension are obtained from 
Ref . [12. 



means that the photons are not able to tunnel to adjacent 
cavities, since too much energy would be needed for this 
process. Hence, the Mott phase can be considered as a 
stable state of frozen light. In superfluid phase, however, 
the polaritons are delocalized on the whole lattice and 
Bose condense in the state of zero momentum. 

We determine the phase boundary of the two- 
dimensional TCL model for zero detuning A = 0. The 
first two Mott lobes with polariton density Up = 1 and 
Up = 2, respectively, are shown in Fig. [2] for distinct val- 
ues of the two- level system number TV = {2, 3, 5, 10}. 
The boundary of the quantum phase transition is given 
by the minimal amount of energy necessary to add (re- 
move) a particle to (from) the system and can therefore 
be evaluated directly from the minimal gap of the single- 
particle spectral function obtained by means of VGA. 
The size of the gap does not depend on the particle 
species the spectral function is evaluate for, since pho- 
tons and atomic excitations are coupled by 5. As already 
mentioned in Section. IIIIl we use the variational param- 
eters X = {ujc, e, i}, which allow for thermodynamic con- 
sistency in the total polariton number. In contrast to 
the results in one dimensior^i^ the lobes are round shaped 
and no reentrance behavior can be observed for increasing 
hopping strength t. For increasing number of two-level 
systems N the width of the Mott lobes with different po- 
lariton density Up becomes more similar, see Tab.]!] The 
critical hopping strength i*, which determines the tip of 
the Mott lobes, depends on both the filling Up and the 
number of two-level systems per cavity N. In particular, 
t*{N) is shrinking for increasing N for a fixed polariton 
density Up. In Fig. [3] we investigate for the first Mott 



lobe, i.e., for rip = 1, the dimensionality dependence of 
the ratio t*{N)/t*{2), which specifies how fast the lobes 
are shrinking with increasing N. To this end, we com- 
pare our VGA results for two dimensions with DMRG 
results for one dimension obtained by D. Rossini et al. 
in Ref. [ij- It can be observed that with increasing A'' the 
lobes are shrinking faster in two dimensions than in one 
dimension. 

The phase boundary at zero hopping can be deter- 
mined analytically, as the model decouples into M single- 
cavity problems, i.e., into M TC systems shifted by 
the chemical potential —fxhp. The TC model has been 
solved exactly for zero detuningiii— and for nonzero 
dctuningjiii^ Since the full analytic solution is involved 
we concentrate here in determining the zero hopping 
phase boundary of the first Mott lobe for zero detuning, 
which is relevant for our data. To this end, we diagonal- 



ize the TC model for 
yields 



0, 1, and 2 polaritons which 



^0 = (10a) 

El = {{lOc -fi)± Vn} (10b) 

E2 = {2{ujc - m) ± v/2(2iV-l), 2iujc - A^)} . (10c) 



The number of cigenstates for the sector of Up polaritons 
is Up + 1 ii Up < N and N + 1 if Up > N,^.L21 The 
phase boundary is evaluated by comparing the ground- 
state energies of adjacent sectors, which yields ^ — lJc ^ 
— yN for the boundary between Up = and Up = 1 and 
/i — cjf, = a/ZV — yj2{2N — 1) for the boundary between 
Up — 1 and Up — 2. This is of course in full agreement 



with our numerical results. In light-matter systems the 
optical nonlinearities arise as the energy which is needed 
to add the first excitation to the system is smaller than 
the one to add the second excitation. This results in a 
repulsive interaction of size 2\/N — y/2{2N — 1), which 
is approximately 1/2^/N for large A^. 



V. EXCITATIONS 

In this section we investigate the excitations of the 
TCL model. In particular, we evaluate single-particle 
spectral functions and densities of states of photons and 
atomic excitations. Furthermore we present the momen- 
tum distribution for both particle species. Based on the 
spectral information we characterize polaritons, which 
are the elementary excitations in light-matter systems. 



A. Spectral properties of photons and atomic 
excitations 

Photon spectral functions ^^''(k, w) evaluated by 
means of VGA for fixed hopping strength t = 0.015 
and zero detuning A = are shown in Fig. 2] for N = 
{2, 3, 5, 10} two-level systems located in each cavity. For 
increasing number of two-level systems N the Mott lobes 
shrink and thus the gap in the spectral function is de- 
creasing. The modified chemical potential /i — Wc is cho- 
sen such that the spectral function is evaluated approx- 
imately in the middle of the Mott lobe. We used the 
variational parameters x = {uc, e, t}, reference systems 
of size L = 4 X 2, and an artificial broadening 0+ = 0.01 
for the numerical evaluation. 

The number of particle and hole bands, present in the 
single-particle spectral function and their approximate 
energies can be already determined from the single-cavity 
solution. For large enough filling {up > N) there are N + 
1 eigenstates in the sectors of rip ± 1 particles which leads 
to A^ -I- 1 particle and hole bands, respectively. However, 
we investigate spectral properties in the first Mott lobe 
(rip = 1) and thus we have to examine the zero- and 
two-particle sectors. The zero-particle sector {up = 0) 
consists of only one state leading to one hole band uJh and 
the two-particle sector (up = 2) consists for all iV > 2 
of three states, which leads to three particle bands w* 
where i G {1, 2, 3}. We choose the order of the bands 
such that the excitation energy increases with increasing 
index i. For clarity Fig. 2] shows only the hole band w/j 
and the lowest-lying particle band w^. The approximate 
location of the particle bands is obtained from the energy 
difference of the eigenenergies of the two particle sector 
and the ground-state energy of the one-particle sector, 
see Eq. (fTOl) , leading to 
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FIG. 4: (Color online) Hole band Uh and the lowest-lying 
particle band ojp of the photon single-particle spectral func- 
tion ^^'''(k, (jj) for fixed hopping strength t = 0.015 and zero 
detuning A = 0. The modified chemical potential and the 
number of two-level systems is (a) ^ — lUc = —1.25, N = 2, 
(b) n-ujc^ -1.6, iV = 3, (c) fi~ujc^ -2.15, N ^ 5, and 
(d) ii-LOc = -3.1, iV = 10. 



ul « (ujc -^l) + V2(2iv-i) + Vn 

Analogously, one obtains for the hole band 



^h 



{ujc -fi)-VN 



w; « (wc - Ai) + viv 



The densities of states of both photons A^^'' (w) as well 
as atomic excitations N'^^{uj) arc shown in the first and 
second row, respectively, of Fig. [5] for identical parame- 
ters as in the case of the single-particle spectral func- 
tion. The left column contains the density of states of a 
small energy window centered around zero, showing the 
low- lying excitation bands uih and cOp. The right column 
contains data for higher excitation energies. The bands 
Wp carry significant spectral weight whereas the bands 
ujp are barely visible for A^ = 2 (w^ w 5.1) and A^ = 3 
{ujp « 6.5). For more than three two-level systems per 
cavity the excitation energy is already larger than the 
maximum energy considered in the plot. The position 
of the bands matches well with the approximate results 
obtained from the single-cavity limit. In the photon den- 
sity of states Np''-{ll!), first row, the spectral weight of the 
low-lying particle band w^ decreases with increasing two- 
level system number A'^, whereas the spectral weight of 
u!p increases. The opposite is true for the atomic excita- 
tion density of states N'^^ (w) , second row. Interestingly, 
due to this behavior photon and atomic excitation densi- 
ties of states become more similar for increasing two-level 
system number A^, which is due to the fact that the en- 
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FIG. 5: (Color online) Density of states of (a) photons N^''{uj) 
and (b) atomic excitations TV^^ (lj) for parameters as in Fig.[4l 
The left column shows contributions from the bands ujh and 
Wp and the right column from the bands with higher energy 
L0„ and some of the to^ . 
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FIG. 6: (Color online) Momentum distribution of (a) photons 
n'' (k) and (b) atomic excitations n°^(k) for parameter as in 
Fig.H 



semble of two- level systems behaves more like free bosons 
for large NM^ 

The momentum distribution of both particle species, 
which can be evaluated with high accuracy by means 
of the Q-matrix formalism^ is shown in Fig. [51 The 
momentum distribution of photons (left panel) and the 
atomic excitations (right panel) do not exhibit major dif- 
ferences. For increasing number of two-level systems N 
and fixed hopping strength t, the tip of the Mott lobe 
is approached, see Fig.[2j From this in turn it follows 
that the density in the center of the Brillouin zone is in- 
creasing for increasing N, which is in accordance with 
our results. 



B. Polaritons 

Here, we investigate properties of polaritons, the ele- 
mentary excitations of the TCL model, which are linear 
combinations of photons and atomic excitations. Our 
goal is to describe these excitations by polaritonic quasi- 
particles added to the A'^p-particle ground state I^/jo). 
Hence, we introduce the polariton creation operators pj^ j, 

for particle excitations and h'^ ^ for hole excitations as 
suitable linear combinations of photons and hard-core 
bosons introduced in Sec. IIIIl 






flla) 



(lib) 



The weights /3" (k) and 7" (k) of the polariton creation 
operators depend on the wave vector k, the quasiparti- 
cle band index a, and the filling Up. The dependence on 
the latter is not explicitly included in the notation as we 
solely focus on the first Mott lobe with particle density 
Up = I. It is important to notice that the hole creation 
operator is neither the adjoint of the particle creation op- 
erator nor its annihilation counterpart. The normalized 
polariton quasiparticle states are generated by applying 
the polariton particle and hole creation operators on the 
A^p-particle ground state. 



IC,k> = 



pIm 1^0) 



Y(V'o|Pa,kpl,klV'o) 



(12a) 
(12b) 



The weights /? and 7 of the linear combination are de- 
termined by maximizing the overlap between the exact 
eigenvectors |^q, ^ ) of the TCL model in the (Np ± 1)- 
particle sector. This yields a generalized eigenvalue prob- 
lem which is used to determine the weights (3 and 7, see 
Ref. [ij for a detailed derivation and discussion. The 
eigenvalue A of the generalized eigenvalue problem spec- 
ifies the quality of the quasiparticle description. More 
specifically, A is bound by the interval [0, 1], where 
A = 1 corresponds to a perfect description by polariton 
quasiparticles, i.e., to maximal overlap between the true 
(iVp± l)-particle states and the polariton states |V'p/;i k)' 
whereas small values of A indicate a modest quasiparti- 
cle description. The generalized eigenvalue problem fixes 
the weights (3 and 7 only upon a constant, which is de- 
termined by the condition that the total spectral weight 
consisting of the spectral weight of photons and atomic 
excitations is conserved. 

For the hole band uth the weights (3 and 7 can be chosen 
freely, as both Ok and 6k applied on the ground state 
with particle density np = I are proportional to the same 
state. Thus we investigate only the weights /3 and 7 for 
the particle bands Wp, Wp, and Wp, which are shown in 
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FIG. 7: (Color online) Photon weights /3 (solid lines) and 
atomic excitation weights 7 (dashed lines) of the polaritonic 
quasiparticle creation operators p]^ ^ for the bands (a) ojp, (b) 
LOp, (c) ujp, and N = {2, 3, 5, 10} two- level systems per cavity. 



Fig.[7]froni top to bottom for N = {2, 3, 5, 10} two-level 
systems per cavity. The photon weight /3 corresponds 
to the solid line and the atomic-excitation weight to the 
dashed line. For the band with lowest excitation energy 
u!p the weights /S and 7 are of opposite sign, whereas the 
sign is equal for the bands uip and ujp. The bands w^ 
and cjp are very well described by the polariton picture 
as A sa 1. Yet, the band with highest energy w^ is very 
poorly represented by the polariton creation operators 
as A « 0.01. For increasing number of two-level systems 
N per cavity the weights of the two constituents become 
similar. This might indicate, as in the case of the spectral 
weight, that the atomic excitations behave for large N 
like bosonic particles. In the case of a single two-level 



system per cavity N = 1, i.e., the Jaynes-Cummings 
lattice model, the asymmetry in the coefhcients is much 
more pronounced than it is here<^ 



VI. CONCLUSIONS 

We presented and discussed the quantum phase tran- 
sition and the excitations of the Tavis-Cummings lattice 
model in two dimensions obtained within the variational 
cluster approach. The Tavis-Cummings lattice model de- 
scribes light-matter systems which contain multiple two- 
level structures in each cavity. Due to this fact the Tavis- 
Cummings lattice model might be easier to realize in the 
experiment than its counterpart, the Jaynes-Cummings 
lattice model, which contains exactly one two-level sys- 
tems per cavity. As a goal for future research, a detailed 
study of cavities with a random number of two-level sys- 
tems might provide further interesting insight into light- 
matter systems. 

In this paper, we determined the quantum phase tran- 
sition delimiting Mott phase, in which polaritons are lo- 
calized in each cavity, from superfluid phase, in which po- 
laritons are delocalized on the whole lattice. We studied 
the dependence of this phase boundary for various two- 
level system numbers per cavity. For increasing number 
of two-level systems the Mott lobes become narrower, 
however, the width of the Mott lobes for distinct filling 
becomes more equal. We also compared the dependence 
of the critical hopping strength, which determines the 
tip of the Mott lobe, on the dimension of the coupled 
cavity system. Additionally to the phase boundary, we 
investigated spectral functions and corresponding densi- 
ties of states. The variational cluster approach allows us 
to extract spectral properties for both photons as well 
as atomic excitations, provided the latter are mapped 
onto hard-core bosons. By investigating the zero-hopping 
limit, which corresponds to investigating a single cavity, 
we determine the number of bands in the spectral func- 
tion and their approximate location. For the first Mott 
lobe there exist three particle bands and one hole band. 
The particle band and the hole band with smallest energy 
are reminiscent of the excitations in the Bose-Hubbard 
model. Particularly, they are also cosine shaped and the 
density distribution of the weight is similar i^i^^"— The 
additional two particle bands lie at considerably higher 
energies. The band with second highest excitation energy 
carries significant spectral weight whereas the one with 
highest energy is barely visible in the spectra. Interest- 
ingly, for increasing two-level system number the weight 
of the photon spectra becomes more and more similar to 
the weight of the atomic excitation spectra. We investi- 
gated the momentum distribution as well, which is rather 
similar for photons and atomic excitations. Yet, for in- 
creasing two-level system number and constant hopping 
strength, a larger density can be observed in the cen- 
ter of the Brillouin zone, which arises due to the fact 
that the Mott lobe is shrinking with increasing number 



of two- level systems. Therefore, for identical hopping 
strength the boundary to superfluid phase is approach- 
ing, which is responsible for the increasing density in the 
center of the Brillouin zone. Finally, we studied the prop- 
erties of polaritons, the elementary excitations in light- 
matter systems. Since we evaluated spectral properties of 
both particle species, we were able to introduce polariton 
quasiparticle and quasihole creation operators as linear 
combinations of photons and atomic excitations. The 
polariton creation operators depend on the wave vector, 
band index and filling. We investigated the photon and 
atomic excitation weights of the linear combination and 
analyzed their dependence on the number of two-level 
systems located in each cavity. 
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